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ABSTRACT 

We measure the luminosity profiles of 16 brightest cluster galaxies (BCGs) at 0.4 < z < 0.8 using high 
resolution F160W NICMOS and F814W WFPC2 HST imaging. The heterogeneous sample is drawn 
from a variety of surveys: seven from clusters in the Einstein Medium Sensitivity Survey (EMSS; Gioia 
& Luppino 1994), five from the Las Campanas Distant Cluster Survey and its northern hemisphere 
precursor (LCDCS; Dalcanton et al. 1997; Gonzalez et al. 2001; Nelson et al. 2001a), and the remaining 
four from traditional optical surveys (Spinrad 1980; Koo 1981; Gunn et al. 1986; Couch et al. 1991). We 
find that the surface brightness profiles of all but three of these BCGs are well described by a standard de 
Vaucouleurs (r^/^) profile out to at least ~ 2re and that the biweight-estimated NICMOS effective radius 
of our high redshift BCGs (r^ = 8.3 ± 1.4 kpc for Ho = 80 km s'^ Mpc-\ = 0.2, D.^ = 0.0) is - 2 
times smaller than that measured for a local BCG sample (Graham et al. 1996). If high redshift BCGs 
are in dynamical equilibrium and satisfy the same scaling relations as low redshift ones, this change in 
size would correspond to a mass growth of a factor of 2 since z ~ 0.5. However, the biweight-estimated 
WFPC2 effective radius of our sample is 18 ± 5.1 kpc, which is fully consistent with the local sample. 
While we can rule out mass accretion rates higher than a factor of 2 in our sample, the discrepancy 
between our NICMOS and WFPC2 results, which after various tests we describe appears to be physical, 
does not yet allow us to place strong constraints on accretion rates below that level. 

Subject headings: galaxies: clusters: general — galaxies: evolution — galaxies: formation — galaxies: 
elliptical 



1. INTRODUCTION 

Brightest cluster galaxies (BCGs) have been studied, 
both photometrically and spectroscopically, in great de- 
tail at low redshift because of their unique characteristics 
and location (e.g., Hoessel et al. 1980; Schombert 1986; 
Postman & Lauer 1995). BCGs are the most massive and 
luminous galaxies in the universe, with central velocity dis- 
persions cr - 300-400 km s'^ (Dressier 1979; Carter 1985; 
Fisher et al. 1995), and typical luminosities ^ 10 L*, where 
L,, = 1.0 X IO^O/i^Lq, (Sandage & Hardy 1973; Hoessel et 
al. 1980; Schombert 1986). Rather than simply being "the 
brightest of the bright," BCGs are drawn from a different 
luminosity function (Dressier 1978) and follow tighter scal- 
ing relations, for example the r^, — fi^ relation (Hoessel et 
al. 1987), than other cluster ellipticals. Their remarkably 
homogeneous luminosity, cry ~ 0.35 mag, (cry ~ 0.30 mag 
when applying corrections based upon environment; e.g., 
Sandage & Hardy 1973; Hoessel et al. 1980), suggests they 
may experience a unique evolutionary history. Because 
BCGs are typically located near the center of the cluster 



potential, they may preferentially accrete material from 
tidally stripped cluster galaxies (Richstone 1976) or merge 
with other cluster members (Ostrikcr & Hausmann 1977). 
Observed correlations between BCG structural parameters 
and luminosity (such as the luminosity-shape relation of 
Hoessel et al. 1980) may be signs of such accretion events 
(Ostriker & Tremaine 1975; Schneider et al. 1983). 

Studies comparing the luminosity of high redshift BCGs 
to their low redshift counterparts find indirect evidence of 
mass accretion in BCGs since z ~ 1 (Aragon-Salamanca 
et al. 1998, hereafter ABK98; Burke et al. 2000; Nelson 
et al. 2001a). The amount of inferred evolution, however, 
is tied to environment. BCGs in clusters with low x-ray 
luminosity (La; < 2.3 x 10^^ ergs s"^) fade less with de- 
creasing redshift than expected from the aging of their 
stellar populations (ABK98; Burke et al. 2000; Nelson et 
al. 2001a). The proposed explanation for this observation 
is that these galaxies are accreting roughly enough mass 
in stars to counter the expected passive evolution - about 
a factor of 2 to 4 increase in mass since z ^ 1 (ABK98). 
On the other hand, the luminosities of high BCGs do 



^ Based on observations made with the NASA/ESA Hubble Space Telescope which is operated by AURA, Inc., under contract with NASA. 

Guest User, Canadian Astronomy Data Center, which is operated by the National Research Council of Canada, Herzberg Institute of 
Astrophysics, Dominion Astrophysical Observatory. 
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match the predictions of passive evolution models (Burke 
et al. 2000; Nelson et al. 2001a), which hmits the possible 
amount of mass these BCGs accrete to less than a factor 
of 2. 

Unfortunately, the inferred mass accretion rates rely 
on the assumption that increases in luminosity are at- 
tributable to increases in mass. A way to test the hy- 
pothesis of accretion in BCGs is to measure the structural 
parameters of the BCGs (mass, size, velocity dispersion, 
surface briglitncss). Locally, BCGs conform to the funda- 
mental plane relations between log Vg— log a — (/u) derived 
by Dressier (1987) for lower luminosity ellipticals (Oegerle 
& Hoessel 1991). Because they lie on the fundamental 
plane, the local BCGs appear to be in dynamical equilib- 
rium and changes in mass should produce corresponding 
changes in size. If we presume that high redshift BCGs 
also satisfy these scaling relations on the grounds that 
high redshift ellipticals arc observed to lie on the funda- 
mental plane (Kelson et al. 1997), then a comparison of 
the sizes of local and distant BCGs should help constrain 
accretion models. Alternatively, highly irregular BCG lu- 
minosity profiles would argue that these galaxies are not 
in dynamical equilibrium. The major obstacle for such 
a study is that measuring structural parameters of high 
redshift BCGs requires high angular resolution imaging. 

We present the first study of the structural parameters 
of high redshift BCGs using Hubble Space Telescope NIC- 
MOS F160W and WFPC2 F814W imaging and search 
for direct signatures of accretion events in these galax- 
ies. The use of infrared data is important because, even 
for our most distant clusters, these wavelengths probe the 
older, more quiescent stellar populations, thereby making 
comparisons to local samples more direct. Without prior 
knowledge, a homogeneous data set that minimizes one's 
sensitivity to potential recent star formation is necessary 
for this study. If the high-redshift BCGs are found to be 
quiescent, relaxed systems, then one might, with caution, 
further explore the WFPC2 archive to enlarge the sample. 
In §2 we describe our heterogeneous sample of 16 BCGs 
at 0.4 < 2 < 0.8 and the HST data. In §3, we describe 
the three different methods we use to measure BCG effec- 
tive radii including GIM2D, a package which we utilize to 
fit the surface brightness profiles of the BCGs. In §4, we 
present the best fitting models to the surface brightness 
profiles and investigate the mass accretion rate of BCGs 
since z ~ 0.5. Finally, we summarize our conclusions in 
§5. 

2. THE DATA 

2.1. HST/NICMOS 

The 16 BCGs we analyze are a heterogeneous sample 
culled from several sources - seven from clusters in the 
Einstein Medium Sensitivity Survey (EMSS; Gioia & Lup- 
pino 1994), five from the Las Campanas Distant Cluster 
Survey and its northern hemisphere precursor (LCDCS; 
Dalcanton et al. 1997; Gonzalez et al. 2001; Nelson et 
al. 2001a), and the remaining four from traditional op- 
tical surveys (Spinrad 1980; Koo 1981; Gunn et al. 1986; 
Couch et al. 1991). The BCGs were observed by the Hub- 
ble Space Telescope with the NICMOS camera and F160W 
filter from June 1997 through March 1998 in the "snap- 
shot" mode as part of program GO 7327 (PI: Dalcanton). 



The images were obtained using the NIC2 camera con- 
figuration (19:'2 X 19:'2 field of view, 0.075" pixel"!) in 
MULTIACCUM mode with NSAMP=16. Observations 
consist of three dithered exposures per object with expo- 
sure times of 256 seconds per dither position for BCGs 
a,t z < 0.6 and 352 seconds per dither position for BCGs 
at z > 0.6. The NIC2 image datasets are listed in Ta- 
ble 1. NICMOS images calibrated "on-the-fly" from the 
HST data archive are generally of inferior quality because 
the automated reduction pipeline, CALNICA. does not 
correct for non-linear bias levels and time- varying changes 
in the dark current level (the so-called "pedestal effect"). 
Therefore, we calibrate the raw science images using the 
NICPIPE routine available in the STSDAS Version 2.2 
IRAF package. This pipeline incorporates the standard 
reduction procedure of CALNICA, but in addition uses 
the tasks BIASEQ and PEDSUB to correct for non-linear 
bias level drift and quadrant-dependent residual bias level. 

We calculate F160W magnitudes (denoted i?i60,AB 
hereafter) in Oke's AB system (Oke & Gunn 1983) ac- 
cording to the following expression: 

Hi6o,AB = -2.5log{PHOTFNU x CR) + 8.9, (1) 

where CR is the count rate and PHOTFNU is the HST 
photometric conversion between countrates and fluxes in 
units of Jy s DN~! in the image headers. Galaxy magni- 
tudes are corrected for extinction using the dust IR emis- 
sion maps of Schlegel et al. (1998) but are not K-corrected. 

The galaxy identified as the BCG was generally not 
known prior to this program. Because the field of view 
of the NICMOS images is small, we do not simply choose 
the brightest galaxy on the image as the BCG. Instead, 
we identify the BCGs on large optical and/or near-IR 
images for all but two clusters MS0302.5-I-1717 and 
MS1333.3-hl725 (Gioia & Luppino 1994). For the LCDCS 
clusters, the BCG is defined as the galaxy with the 
brightest total /-band magnitude located within 350h~^ 
kpc of the cluster center. To reduce contamination by 
bright blue foreground galaxies, we exclude BCG candi- 
dates whose colors are 0.4 mag bluer than the location 
of the red envelope in the cluster color-magnitude dia- 
gram (sec Nelson et al. 2001a for details). 1100+4620 
does not have /-band imaging and therefore the BCG is 
determined using ground-based Kp imaging. We identify 
the BCGs in CL0016-I-16, J1888.16CL, CL0317-hl521, and 
CL1322-)-3027 using the cluster contour plots and /T-band 
photometry files published by Aragon-Salamanca et al. 
(1993). Five of the seven EMSS BCGs are selected us- 
ing various optical images and identifications published 
in the literature (Luppino & Gioia 1992; Clowe 1998; 
Luppino & Gioia 1999; Clowe et al. 2000). Finally, we 
could not find BCG identifications in the literature for 
the remaining two EMSS clusters (MS0302.5-I-1717 and 
MS1333.3-M725). The NICMOS images for both of these 
clusters contain a large, dominant galaxy in the field of 
view which we define to be the BCG. We caution, how- 
ever, that the BCG identifications for these two clusters 
are less certain. 

2.2. HST/WFPC2 

We also make use of archival HST/WFPC2 images taken 
in the F814W filter available for 9 BCG's in our NIC- 
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MOS sample. The datasets are listed with their exposure 
times in Table 1. These archival images were recalibrated 
"on-the-fly" through the Canadian Astronomy Data Cen- 
ter (CADC) standard pipeline. WFPC2 images covers a 
larger field of view (36!'8 x 36 "8 for the Planetary Camera 
and 80" x 80" for each of the three Wide-Field Camera de- 
tectors) than NIC2. kt z = 0.6, the F814W filter samples 
the rest-frame wavelength range 4578-5676 A. As a com- 
parison, the NIC2 F160W filter at that redshift samples 
the rest-frame wavelength range 8750-11250 A. 

3. ANALYSIS 

One of the main problems in measuring the sizes of high- 
redshift brightest cluster galaxies are the systematic errors 
in the measurement of the sky. This problem is exacer- 
bated by the fact that some BCGs locally exhibit steep 
(Sersic index n > A) profiles and a significant amount of 
flux in their broad outer wings (Graham ct al. 1996). Sky 
errors are of particular concern for the NIC2 detector be- 
cause it covers a relatively small sky area, which can lead 
to contamination of flat-fielding and sky subtraction. Be- 
cause of the uncertainties involved in the sky determina- 
tion, we decided to measure BCG sizes in three differ- 
ent ways: curve of growth with circular apertures (Sec- 
tion 3.1), Point- Spread-Function (PSF)-convolved surface 
brightness models (Section 3.2), and isophotal ellipse fit- 
ting (Section 3.3). 

3.1. Curve of Growth 

For each BCG in our sample, we measure NIC2 and 
WFPC2 circular half-light radii by computing the galaxy 
flux within concentric apertures and interpolating this 
curve of growth to flnd the radius that contains half of 
the total flux given by the asymptotic flux value at large 
radii. The curves of growth were measured directly from 
the images and not from surface brightness models such as 
the ones obtained in Section 3.2. Typical maximum radii 
for the NIC2 and WFPC2 curves of growth are 6.8-9.0 
arcseconds (90—120 pixels) and 20 arcseconds (200 pixels) 
respectively. The sky background is manually adjusted un- 
til the integrated flux is asymptotically fiat over a range 
of at least 30 pixels. Flux from neighboring galaxies is 
carefully excluded from the curves of growth using seg- 
mentation images produced with the SExtractor galaxy 
photometry package (Bertin & Arnouts 1996). The actual 
BCG flux in each excluded pixel, which would have been 
included in the curve of growth in the absence of neigh- 
bors, is calculated using its symmetric counterpart about 
the center of the BCG. The NIC2 curve-of-growth BCG 
half-light radii (in arcseconds and in kpc) is presented in 
Table 2 (Columns 3 and 4). We list three values for each 
BCG because a different curve of growth was obtained in- 
dependently for each of the three N1C2 dither positions. 
The typical variations in half-light radii from one dither 
position to another are - 0''l-0"2. The WFPC2 curve- 
of-growth half-light radii are presented in Table 3. 

3.2. BCG Surface Brightness Profile Fitting 

The BCG surface brightness proflles are flt using 
the GIM2D package, which is designed to perform 2D 
disk-|-bulge deconvolutions of low signal-to-noise images 
of distant galaxies (Simard 1998; Simard et al. 2001). 



The GIM2D bulge component proflle is a Sersic proflle 
(Sersic 1968) of the form: 

S(r) = Eeexp{-6[(r/re)i/" - 1]}, (2) 

where E(r) is the surface brightness along the semimajor 
axis at radius r, r,, is the semi-major axis effective radius, 
and Se is the effective surface brightness. The parameter 
b is set equal to 1.9992n — 0.3271 to ensure that re is the 
projected radius enclosing half of the light. A standard de 
Vaucouleurs profile is obtained by setting n = 4. The fiux 
interior to the radius r is given by 

F{r)=2nnE,rl^j{2n,b), (3) 

where 7(2n,.T) is the incomplete Gamma function with 
X = h{r /r^)^/'^ . Equations 2 and 3 are given in their circu- 
larly symmetric form for the sake of simplicity, but GIM2D 
does consider the intrinsic (i.e. before PSF convolution) 
ellipticity of the bulge by including bulge ellipticity e and 
position angle (f)b as additional parameters. 

The GIM2D disk component is a simple exponential pro- 
file of the form: 

S(r) = T.Qexp{-r/rd). (4) 

Eg is the face-on central surface brightness, and is the 
semi-major axis exponential disk scale length. The disk is 
assumed to be infinitely thin. The total flux in the disk is 
given by: 

Fdisk - 2^r2So. (5) 

The projected surface brightness distribution of the disk 
inclined at any angle i is calculated by integrating Equa- 
tion 4 over the area in the face-on disk plane seen by each 
projected pixel. The position angle (pd of the disk is also 
a fltting parameter. 

Mosaicing dithered images with non-integer pixel shifts 
requires flux interpolation between pixels and corrections 
for geometric distortion. For undersampled HST images 
such as ours, these processes severely distort the shape of 
the PSF and destroy the uncorrelated noise characteris- 
tics of theoriginal images. To circumvent these potential 
difficulties, GIM2D allows for fitting multiple, dithered in- 
dividual images simultaneously. We use this mode for the 
NIC2 data. The WFPC2 images are also undersampled, 
but because our archival images have integer pixel shifts 
they are combined without information loss, and the com- 
bined image is analyzed with GIM2D. 

HST point-spread- functions vary spatially across the im- 
age, and additionally NICMOS PSFs vary temporally. A 
PSF must therefore be individually generated for each ob- 
ject of interest. We use the software package TINYTIM 
V5.0c (Krist 1993) to generate PSFs that include both 
types of variation for each BCG position and observation 
date (for NIC2). Our NIC2 PSFs are i;'9 x 1'.'9, and 
our WFPC2 PSFs are 2:'4 x 2 ('4. Both types of PSFs 
are oversampled by a factor of 5. To compare directly to 
undersampled HST data, GIM2D creates a surface bright- 
ness model on an oversampled grid, convolves it with the 
corresponding oversampled PSF, shifts it to a subpixel po- 
sition specified by the values of the fitting parameter dx 
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and dy and robins the resulting profile to the detector's 
resolution. 

We use SExtractor segmentation images to isolate the 
BCGs from their neighbors. As SExtractor performs 
source detection, it deblcnds sourc;c;s using flux multi- 
thresholding (Bertin & Arnouts 1996). The result of de- 
blending are segmentation images in which pixels belong- 
ing to the same object are all assigned the same flag value, 
and sky background pixels are flagged by zeroes. GIM2D 
fits arc performed on all pixels flagged as the target ob- 
ject or background in the SExtractor segmentation im- 
ages. The rationale behind not just fitting "object" pixels 
is that considering only the pixels interior to the bound- 
ary between object and sky, which is sharply delineated 
by the isophote corresponding to the detection threshold, 
would exclude information contained in the pixels below 
that threshold. 

We fit the NIC2 BCG images with pure de Vaucouleurs 
(n = 4) bulge and full de Vaucouleurs bulge + exponential 
disk models. Although the sky level can be a free parame- 
ter in GIM2D, we fix the sky level to the value determined 
from the curve of growth analysis of Section 3.1 because of 
the relatively large angular extent of the galaxies and the 
non- uniformity of the background. The results from the 
GIM2D NIC2 fits are shown in Table 2. The GIM2D half- 
light radii are measured along the semi-major axis (SMA) 
of the BCG. SMA radii can be converted to "circular" radii 
according to rdrc — fsma -^/l — Cbcg whore e^cg is the ellip- 
ticity of the BCG. We list both SMA and circular radii 
in Table 2. Residual images are visually inspected to de- 
termine which model (pure bulge or bulge+disk) yields a 
better fit. The better model is listed in Column 13 ("PB" 
= Pure Bulge, "BD" = Bulge-hDisk and "EQ" = Equally 
Good) of Table 2. Although seven BCGs are fit better by 
a bulge-|-disk model, this result does not imply that actual 
disks are present in those BCGs. The residual images sug- 
gest that an intrinsic change in ellipticity (i.e. over what 
one would expect from pure PSF smoothing) and position 
angle as a function of radius may be responsible for the 
preference of a bulge-l-disk model. 

We fit the deeper and larger WFPC2 images with three 
different models: a pure de Vaucouleurs bulge, a de Vau- 
couleurs bulge + exponential disk and a pure Sersic of pa- 
rameter n. We allow n to vary between 1 and 6. Two sets 
of fits are performed: one with the sky fixed to the value 
obtained from the curve of growth analysis and one with 
the sky left as a free parameter. The results are shown 
in Table 3. There are two lines per BCG in Table 3. The 
first and second lines give the results of the fixed and fioat- 
ing background fits respectively. The half-light radius of 
the best model, as determined through visual inspection 
of the residual images, is underlined for each BCG. The 
advantage of using so many different models is the ability 
to sec the systematic effects on the measurement of the 
BCG half-light radii introduced by the choice of a particu- 
lar model. The significant disagreements between models 
arises because the half-light radius depends on the mea- 
sured total flux, which depends on the details of the poorly 
constrained outer regions of the galaxy. Table 3 lists six 
radii for each BCG, and the variations among these six 
values reflect errors due to the modeling of the sky and 
the flux in the outer wings of the galaxy. 



3.3. Elliptical Isophote Fitting 

Elliptical isophote fitting provides yet another way to 
measure BCG sizes and to demonstrate how bulge+disk 
model fits can introduce a disk component to model a 
pure elliptical galaxy with ellipticity and position an- 
gle radial variations. We use the task ELLIPSE in the 
IRAF/STSDAS/ISOPHOTE package to fit isophotal el- 
lipses to the WFPC2 images of seven out of nine BCGs. 
We could not acceptably fit MS1621 and MS0451. The 
"BCG" in MS1621 is actually a close pair of early-type 
looking galaxies, and the MS0451 BCG is adjacent to an- 
other large galaxy. 

In this fitting, we allow the isophote center, ellipticity, 
and position angle to vary with radius. In contrast, in- 
troducing radial variations is not practical in model fit- 
ting because of the large number of degrees of freedom. 
For example, the GIM2D model has one center, one el- 
lipticity and one position angle per component (bulge or 
disk). The isophotal ellipse fitting results are presented in 
Table 4. The Table presents the difference in ellipticity 
between the outer and inner isophotes, Ae(ISO); the dif- 
ference in ellipticity between the outer (disk) component 
and the inner (bulge) component, Ae(GIM2D), the dif- 
ference in isophote position angle between the outer and 
inner isophotes, A(/)(ISO), the position angle difference be- 
tween the disk and bulge components, A(/!)(GIM2D), and 
the pixel shifts in the isophote center between the outer 
and inner isophotes. Ax and Ay. The close correlations 
between Ae(ISO) and Ae(GIM2D) and between A(^(ISO) 
and A(/)(GIM2D) suggest that bulge-|-disk models fit some 
of the BCG profiles better than bulge-only models because 
the extra freedom in the two component models is being 
used to fit the radial variations in the ellipticity and posi- 
tion angle. 

3.4. NIC2 Curve- of- Growth Simulations 

To test how well we recover true half-light radii from 
NIC2 curves of growth, we perform a set of simulations. 
We choose five BCGs with both NIC2 and WFPC2 images: 
MS1137, CL1322, MS2053 and CL0016. We create three 
types of NIC2 simulations for each BCG. For the first sim- 
ulation, we take the NIC2 best-fit pure r^/^ GIM2D model 
image, and we insert it in another location on the NIC2 
image. This type of simulations is labeled as "NIC2 r^^'*" 
in Table 5. For the second simulation (labeled as "NIC2 
B-hD" in Table 5), we take the NIC2 best-fit bulge-|-disk 
GIM2D model, and we insert it in another location on 
the NIC2. For the third simulation (labeled as "WFPC2 
7.1/4" Table 5), we use the half-light radius of the best-fit 
WFPC2 model to create a pure r^/^ model. We normalize 
this model to have the measured fiux, and we insert it in 
the NIC2 image. We analyze all three types of simulations 
using the same curve-of-growth procedure that we use to 
measure circular half-light radii from the real data. See 
Table 5 for the results. No model fits are involved. We 
compare the radii recovered from the curve-of-growth to 
the circular input values. Fr is the ratio of the recovered 
and input radii, and Ff is the ratio of the recovered and 
input total fiuxes. Using the results from all three types 
of simulations, we find that the NIC2 half-light radius is 
systematically underestimated by a mean of 20%. The 
cause of the problem is not simply the lack of sky in the 



Constraints on Size Evolution of Brightest Cluster Galaxies 



5 



NIC2 images because one of the galaxies with the largest 
WFPC2 half-light radii (MS1137) is recovered well. 

4. DISCUSSION 

4.1. Are High-Redshift BCGs in Dynamical Equilibrium? 

The presence of significant distortions indicating recent 
major mergers and/or an ellipticity distribution that is sig- 
nificantly different from the local BCG population would 
be clear morphological signs that high-redshift BCGs are 
not in dynamical equilibrium. However, the surface bright- 
ness profiles of our 16 BCGs arc quite regular (Figure 1). 
To generate the plotted profiles, we fit elliptical isophotes 
to the data using the IRAF task ELLIPSE. We use the 
best fitting model galaxy parameters listed in Table 1 
to compute the slope and zero-point of the correspond- 
ing PSF- deconvolved de Vaucouleurs profile. The semi- 
major axis radii are converted from arcsec to kpc assum- 
ing Ho = 80 km s-i Mpc-\ fi^ = 0.2, and 9.^ = 0.0. 
The BCG surface brightnesses have been corrected for 
cosmological (1 + z)^ dimming, but arc not K-corrected 
(K-corrcctions for an elliptical galaxy in the TJ-band are 
typically ~0.2-0.4 mag for this redshift range). Figure 1 
demonstrates that all but three of our BCGs are well de- 
scribed by a de Vaucouleurs profile beyond the resolution 
limit and out to at least ~2 NIC2 Ve- The notable excep- 
tions are MS0302.5-M717, CL0016, and 1041-^4626. The 
BCG MS0302.5-M717 is actually best-fit by a pure expo- 
nential disk profile. We note, however, that this BCG is 
selected directly from the NICMOS images and therefore 
its identification is less certain. The BCG CL0016 has 
a prominent dust lane, and the BCG 1041-1-4626 lies at 
the corner of a frame where the local background varies 
significantly. Additionally, we search for azimuthal asym- 
metry in the BCG profiles by subtracting the best fitting 
2D galaxy model from the data. With the exception of 
the BCG in MS0302.5-I-1717, which is not well fit by an 
profile, and the BCG in CL0016, which shows a dust 
lane, we find that the BCG images do not have significant 
2-dimensional residuals (the best fitting models typically 
yield ~ 1). Consistent with this finding, our high red- 
shift BCGs have similar ellipticities (0.0 < e < 0.5) to 
their local counterparts (Graham ct al. 1996). We con- 
clude that the BCGs exhibit no signs of being far from 
dynamical equilibrium. 

4.2. Comparison with Local BCG Sizes 

The good agreement between our BCG surface bright- 
ness profiles and a standard de Vaucouleurs profile in con- 
junction with the lack of any significant 2-dimensional 
residuals argues against recent major accretion events. 
However, cumulative minor merger events over time may 
manifest themselves not by gross morphological distor- 
tions, but by the gradual enlargement of the central galaxy. 
The simplest way to quantify the sizes of these galaxies is 
through Tg. In Figure 2 we show the distribution of semi- 
major axis Te from the GIM2D pure de Vaucouleurs bulge 
fits to our BCGs {black histogram). To avoid "double- 
weighting" MS1621 because it is a binary BCG, we only 
include BCG MS1621a. The choice of BCG in MS1621 is 
inconsequential for this discussion because both galaxies of 
the "BCG pair" have similar effective radii. We find that 
our high redshift BCGs have a biweight-estimated effective 



radius of 8.3 ±1.4 kpc in the NICMOS images. We com- 
pare to the local sample of Graham et al. (1996; shaded 
histogram) who fit the surface brightness profiles of BCGs 
from 119 Abell clusters at z < 0.05 in the i?-band with 
a standard de Vaucouleurs profile (i.e. n — A) assuming 
Ho = 80 km s-i Mpc"!. Their low redshift BCGs have 
a biweight-estimated effective radius (measured along the 
semi-major axis as for our sample) of 16.1 ± 1.7 kpc which 
is ~ 2 times larger than that of our NIC2 high redshift 
sample. To quantitatively assess the significance of this 
difference in r^, we randomly select subsets of 16 BCGs 
from the Graham et al. (1996) local BCG sample and com- 
pute a biweight-estimated Tg for each subset. Using 500 
Monte Carlo realizations of the high redshift BCG sample 
drawn from the low redshift sample, we find that < 3% 
(13 of 500) of the realizations have a biweight-estimated 
re less than 8.3 kpc. However, the result is qualitatively 
different if we use the F814VF re's of the 9 BCGs with 
WFPC2 images. The biweight-estimated r^ for that sam- 
ple is 18.2 ± 5.1 kpc — larger than, and fully consistent 
with, the radii of the local sample. 



4.3. Are Optical and Infrared Radii Expected to be the 

Same? 

The discrepancy between the results for the NIC2 and 

WFPC2 radii of our high-redshift BCG sample is either 
due to systematic errors in the NIC2 measurements or to 
a real feature of BCGs. If systematic errors are to blame, 
they would have to be substantially larger than the simu- 
lations in Section 3.4 would suggest. Is it possible that ef- 
fective radii in the rest-frame range 4578—5676 A (F814W 
ai z = 0.6) are systematically different from effective radii 
measured in the rest-frame range 8750-11250 A (F160W 
&t z = 0.6)? Unfortunately, there are no large local sam- 
ples of BCGs with both optical and near-infrared photo- 
metric parameters. However, a tantalizing hint may come 
from the E/SO sample of Pahre (1999). Pahre measured 
global photometric parameters for a sample of 341 nearby 
early-type galaxies in the near-infrared ii'-band. The ma- 
jority (85%) of the galaxies in the Pahre sample reside in 
13 rich clusters with additional galaxies drawn from loose 
groups (12%) and the field (3%). The JsT-band sample was 
complemented with optical F-band data from the litera- 
ture. In Figure 3 we plot the T^-band effective radius versus 
the K-hand effective radius for the 273 galaxies with both 
measurements. The dashed line is the one-to-one line. It 
is immediately obvious that the l^-band radii are system- 
atically larger than the if-band ones, and the discrepancy 
increases with increasing radius. At log = 1.5 in K, the 
V radii are about 1.8 times larger than the K radii. As 
explained in Pahre et al. (1998), this difference is due to 
the presence of color gradients. Given that BCGs are a 
distinct class of objects, it is unclear whether BCG radii 
coiild exhibit the same type of wavelength dependence as 
the E/SO's. However, if they do, then the discrepancy be- 
tween our NIC2 and WFPC2 radii would be explained, 
and the actual size evolution of the BCGs between 2 = 
and 0.6 would be bounded by the WFPC2 images, which 
are bluer in the rest-frame than the local sample, and the 
NICMOS images, which are redder in the rest-frame than 
the local sample. 
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Fig. 1. — The surface brightness profiles of our 16 BCGs ordered by increasing redshift. The filled circles are the data and the line is the 
best fitting de Vaucouleurs (r^/*) profile. The vertical mark in the upper left corner denotes the angular resolution limit (i.e. FWHM of the 
point spread function) of the NICMOS images. All but three of our BCGs are well described by a de Vaucouleurs profile beyond the resolution 
Umit and out to at least ~2 NIC2 r^- The notable exceptions are MS0302.5+1717, CL0016, and 1041+4626. The BCG MS0302.5+1717 is 
actually best-fit by a pure exponential disk profile, the BCG CL0016 has a prominent dust lane, and the BCG 1041+4626 lies at the corner 
of a frame where the local background varies significantly. 




50 100 150 

(kpc) 



Fig. 2. — The distribution of NIC2 semi-major axis eflfective radius, re, from the GIM2D pure de Vaucouleurs bulge fits to our high-redshift 
BCG sample {black histogram) compared to the local Abell sample of Graham et al. (1996) {shaded histogram). 
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Fig. 3. — Log V eflfective radius and log K eflfective radius in arcseconds for 273 E/SO galaxies from the sample of Pahre (1999). The dashed 
line is the one-to-one line. 



4.4. Accretion Rates 

Observations of BCGs show that the relation between 
log a and M (i.e., the Faber-Jackson relation; 1976) is 
rather flat. i.e. BCGs are ovcrhiminous for their velocity 
dispersions (Malumuth & Kirshner 1981, 1985; Oegerle 
& Hoessel 1991). This result has been interpreted as 
support for dissipationless merger scenarios for BCGs in 
which their masses, luminosities, and effective radii in- 
crease, but their velocity dispersions remain constant be- 
cause oc M/R. Because our high redshift BCGs do 
not exhibit any signs of recent major interactions, we pre- 
sume that they are in dynamical equilibrium and that they 
follow the same scaling relations as local BCGs. Conse- 
quently, if their mass scales directly with r,. , wc would infer 
from the comparison of NIC2 and local effective radii that 
BCGs as a class have doubled in mass since z ~ 0.5 — 0.6. 
Using the semi-analytic model predictions of both the 
Durham (Cole et al. 1994; Baugh et al. 1996) and Mu- 
nich (Kauffman et al. 1993) groups, ABK98 present the 
increase in BCG mass as a function of redshift. For CDM 
models with il^ = 1 our inferred mass accretion rate since 
z 0.5 is in excellent agreement with the model predic- 
tions. For low density CDM models with our adopted 
cosmology (fi,„ = 0.2), they predict slightly less mass ac- 
cretion by BCGs since z ~ 0.5 (a factor of ~1.5 increase 
in mass), but given the large observational uncertainties 
our findings arc still in agreement. 

On the other hand, the effective radii from the WFPC2 
images show no significant difference from the local BCG 
sample and would thus be consistent with very little or no 
mass accretion since 2; = 0.5 — 0.6. If this is indeed the 
case, then wc may be faced with the continuing puzzle of a 
class of galaxies with sizes and magnitudes that show little 
or no discernible evolution with redshift even though their 
stellar population must at least be passively evolving (and 



are observed to be passively evolving; Nelson et al. 2001b). 
Part of the possible resolution to this puzzle lies in the 

fact that evolution is tied to environment. The majority 
of our BCGs reside in clusters that fall in the high X-ray 
luminosity category {L^ > 2.3 x 10^** ergs s~^). As noted 
in the introduction, the luminosities of high Lx BCGs do 
match the predictions of passive evolution models (Burke 
et al. 2000; Nelson et al. 2001b), which limits the possible 
amount of mass these BCGs accrete to less than a factor 
of 2. Despite their inherent systematic uncertainties, our 
NIC2/WFPC2 results are sufficient to rule out accretion 
rates higher than a factor of 2 in high BCGs, but our 
results cannot be used to place constraints on accretion 
rates below that level. 

We emphasize that our inferred BCG mass accretion 
rates do not rest on the assumption that the same galaxy 
identified as the BCG at high redshift remains the BCG 
to the present day. Indeed the brightest galaxy in any 
given cluster may very well change over time, but this 
process is part of the evolution one aims to study. The 
value of BCGs in evolutionary studies is their relatively 
unambiguous selection criteria as a class in both the data 
and simulations. 

5. SUMMARY 

Using high angular resolution F160W NICMOS imag- 
ing, we measure the luminosity profiles of 16 BCGs at 
0.4 < z < 0.8. The heterogeneous sample is drawn from 
existing optical and x-ray surveys published in the liter- 
ature (Spinrad 1980; Koo 1981; Gunn et al. 1986; Couch 
et al. 1991; Gioia & Luppino 1994; Gonzalez et al. 2001). 
Archival WFPC2 F814W images were also available for 9 
of the 16 BCGs. Our primary conclusions are that: 
1) the NICMOS surface brightness profiles of high redshift 
BCGs are well described by a de Vaucouleurs profile (r^/^) 
out to a radius of at least ~ 2 NIC2 rg, 
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2) the NICMOS effective radii of our high redshift BCGs 
are, on average, ~ 2 times smaller than those of local 
samples (Graham et al. 1996) for Hq = 80 km s~^ Mpc~^, 
flm = 0.2, and fl\ = 0.0. Using simple dynamical equi- 
librium arguments, this result suggests that BCGs have 
increased in mass by a factor of ~ 2 since z ^ 0.5. 

3) the WFPC2 effective radii are fully consistent with 
those of the local sample, which would suggest little or no 
mass accretion. The discrepancy between the NICMOS 
and WFPC2 results is cither due to systematic errors in 
the NICMOS measurements or to a real physical property 
of BCGs such as color gradients. If systematic errors are 
to blame, they would have to be substantially larger than 
the 20% systematic radius error suggested by our sim- 
ulations. Supporting our suggestion that this difference 
is not due to a systematic error, we note that the optical 
effective radii of local E/SO galaxies are larger than their 
near-infrared values due to color gradients. If the effective 
radii of BCGs exhibit a similar wavelength dependence, it 
would reconcile our NICMOS and WFPC2 radii. Given 
that our WFPC2 and NICMOS data are both bluer and 
redder in the rest-frame than the R-band local sample, the 
results from the two data sets bound the possible evolution 
of Ve, and hence the accretion. 

Although this sample marks a significant improvement 
in the quality of data available for high redshift BCGs, our 
sample is too small to fully resolve systematic uncertainties 
and to confidently constrain mass acc;retion rates below a 
factor of 2. However, our work demonstrates the feasi- 



bility of using high resolution imaging to measure BCG 
structural parameters at z > 0.5. With a sufficiently large 
sample, complemented with some velocity dispersion mea- 
surements to determine whether the scaling relations hold, 
such an approach has the potential to definitively address 
the issue of mass accretion by BCGs since z ~ 1. 
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Table 1 



HST NICMOS AND WFPC2 BCG Image Datasets 



Cluster 


a 


<5 


Total Integration 


CADC Association 


ID 


(J2000.0) 


(J2000.0) 


Time (s) 


Dataset ID 








F160W 


F814W 


F160W 


F814W 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


CL0016+16 


00:18:33.64 


+16:26:15.13 
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8400 


N42B38010 


U2C40101A 


J1888.16CL 


00:56:56.78 


-27:40:31.26 


768 


8400 


N42B39010 


U2C40501A 


MS0302. 5+1717 


03:05:18.19 


+ 17:28:23.35 


768 




N42B48010 




MS0302. 7+1658 


03:05:31.71 


+ 17:10:02.64 


768 


2400 


N42B49010 
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CL0317+1521^ 


03:20:01.52 


+15:31:57.65 


256 


2500 


N42B40010 


U58J0101A 


MS045 1.6-0305 


04:54:10.83 


-03:00:56.87 


768 


10400 


N42B50010 


U3060101A 


1041+4626 


10:41:03.81 


+46:26:36.24 


1056 




N42B05010 




1100+4620 


11:00:57.40 


+46:20:37.88 


768 




N42B15010 




1139-1217 
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768 
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MS1137.5+6625 


11:40:23.86 


+66:08:18.47 


1056 


14400 


N42B52010 


U3060301A 


1147-1252 


11:47:17.30 


-12:52:09.00 


768 




N42B27010 




1230+4621*^ 


12:30:16.26 


+46:21:17.27 


512 




N42B07010 




CL1322+3027 


13:24:48.83 


+30:11:38.16 


1056 


14000 


N42B43010 


U2840501A 


MS1333.3+1725 


13:35:47.14 
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768 
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N42B56010 


U4F30601A 



"^Two of the three dithered exposures for this object had pointing errors 

''One of the three dithered exposures for this object was lost due to data corruption. 
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Table 2 

NIC2/F160W BCG Size Measurements 
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'^One of the three dithered exposures for this object was lost due to data corruption. 
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Table 3 

WFPC2/F814W BCG Size Measurements 
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4.26 


22, 


.06 


3, 
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6, 
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34, 
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.98 
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4, 
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4.77 
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6, 
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,2 
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1, 
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7. 
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1, 
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7, 


,43 


3 


,4 
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8, 
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1, 


.93 


8, 


.70 
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7, 


.84 


1.23 


5, 


.54 


1, 


.05 


4. 


,73 


1 


.02 


4, 


,60 


2, 


,5 


MS2053 


1, 


.77 


9.33 


2.32 


12, 


.23 


2, 


.05 


10, 


.81 


2.33 


12, 


.28 


2.02 


10, 


.65 


2, 


,19 


11, 


,54 


1 
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10, 
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3. 


,8 
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12 
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2, 
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10, 


.54 
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8, 


.01 


1.32 


6, 
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1, 


,93 


10, 


,17 


1 


.69 


8, 
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3, 
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Table 4 



WFPC2 BCG IsoPHOTAL Ellipse Fitting 



BCG 


Fitting Range 


re 


Ae 


Ae 


Acj) 


A4> 


Ax 


Ay 




in r^/^ 


sma 


(ISO) 


(GIM2D) 


(ISO) 


(GIM2D) 


pixels 


pixels 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


CL0016 


1.1 - 2.8 


8.35 


+0.14 


+0.19 


-23 


-31 


+3 


-2 


J1888 


1.0 - 2.1 


5.71 


+0.10 


0.00 


-25 


-30 








MS0302.7 


1.1 - 3.0 


7.02 


+0.15 


+0.39 


+20 


+10 





-3 


CL0317 


1.5 - 2.3 


0.28 


+0.10 


+0.25 





+ 2 





+2 


MS1137 


1.1 - 2.6 


6.32 


+0.14 


+0.26 


-30 


-47 








CL1322 


1.1 - 2.3 


1.31 


-0.05 


0.00 


-20 


+20 


+2 


+2 


MS2053 


1.3 - 2.6 


2.35 


-0.10 


-0.20 


+24 


+60 
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Table 5 

NIC2 BCG Size Measurement Simulations 



BCG 


Model 


Input Thi 


Recovered 




Ff 




Type 


sma 


circ 


COG ru 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


CL0016 


NIC2 r^l^ 


2.89 


2.26 


0.95 


0.42 


0.51 




NIC2 B+D 


1.28 


0.92 


1.02 


1.10 


0.89 




WFPC2 r^l^ 


4.86 


3.80 


1.12 


0.30 


0.38 


MS0302.7 


NIC2 


2.38 


2.26 


1.73 


0.77 


0.82 




NIC2 B+D 


1.61 


1.42 


1.72 


1.21 


1.07 




WFPC2 r^l^ 


3.46 


3.28 


3.97 


1.21 


1.16 


MS1137 


NIC2 


1.89 


1.81 


1.53 


0.85 


0.88 




NIC2 B+D 


1.92 


1.86 


1.56 


0.84 


0.89 




WPFC2 


3.72 


3.55 


2.90 


0.82 


0.95 


CL1322 


NIC2 


0.78 


0.57 


0.46 


0.81 


0.66 




NIC2 B+D 


0.76 


0.55 


0.42 


0.77 


0.63 




WPFC2 


2.70 


1.97 


0.61 


0.31 


0.28 


MS2053 


NIC2 


1.82 


1.62 


1.42 


0.88 


0.98 




NIC2 B+D 


1.77 


1.70 


1.61 


0.94 


1.04 




WFPC2 r^l^ 


2.52 


2.24 


1.68 


0.75 


0.90 



